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^. ' Abstract 



We give an algorithm for obtaining expansions of massive two-loop Feynman 
graphs in powers of the external momentum around a finite, nonzero value of 
C^ ' the momentum. This is based on our general two-loop formalism to reduce 

massive two-loop graphs with renormalizable interactions into a standard set 
of special functions. After the algebraic reduction, the final results are obtained 
by numerical integration. We apply the expansion algorithm to treat the top- 
dependent corrections of 0{g^as) to the b quark self-energy and extract its 
momentum expansion on-shell. 
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Abstract 

We give an algorithm for obtaining expansions of massive two-loop Feynman 
graphs in powers of the external momentum around a finite, nonzero value of 
the momentum. This is based on our general two-loop formalism to reduce 
massive two-loop graphs with renormalizable interactions into a standard set 
of special functions. After the algebraic reduction, the final results are obtained 
by numerical integration. We apply the expansion algorithm to treat the top- 
dependent corrections of O^g^ag) to the h quark self-energy and extract its 
momentum expansion on-shell. 

Analyses of electroweak data based on two-loop radiative corrections have become 
necessary due to the precision level attained already. The experimental precision is 
expected to reach even deeper at the next generation of colliders in the near future. 
This motivated in recent years an enormous amount of work in handling massive 
two-loop Feynman graphs. In its general form, this is a notoriously difficult problem. 
The essence of the difficulties is that these graphs usually lead to unknown special 
functions. 

In a previous paper [|| we gave a general formalism for evaluating two-loop 
Feynman diagrams for arbitrary kinematical variables, based on analytical-numerical 
methods. By using this framework, the internal mass structure and the external 
kinematics of the physical process under consideration are fully respected. Radiative 
corrections are obtained at the end by high accuracy numerical integration. 

In this letter we would like to discuss the use of our general massive two-loop 
framework of ref. |1[ for obtaining external momentum expansions around a finite 
value. There are situations where such momentum expansions around a finite value 
are needed in the evaluation of radiative corrections, most notably for mass and wave 
function renormalization constants. However, we note that the algorithm described 
in this letter can produce the full momentum expansion around a finite value; it works 
systematically for terms of higher order in the momentum as well. 
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We refer to the important process Z ^ bb |0|-|T^ to illustrate in some detail the 
technical aspects. A piece of the analysis of this process involves the calculation of 
the b quark wave function renormalization on-shell. This is a good example for the 
general problem of expanding an amplitude at a certain kinematical point, and we 
would like to use this short communication to treat this problem. Let us begin with 
a broader discussion. 

In a perturbative calculation, one is faced with a set of momentum space integrals, 
where the external momenta appear at various parts of the integrands. The depen- 
dence on the invariants {s} constructed from these external momenta is implicit at 
this stage. One may introduce Feynman parameters and integrate over the internal 
momenta, the result of which will be some Feynman parameter integrals, where the 
dependence on {s} is made explicit. In principle, one can then expand the integrands 
around some fixed values {sq}. 

At the one-loop level, this procedure can be easily followed, since the one-loop 
scalar integrals are relatively simple. At higher loops, there are at least two com- 
plications to follow this route. First of all, due to the proliferation of diagrams, one 
must have a systematic procedure. The Feynman parameter integrals due to different 
diagrams are most likely quite different, and for avoiding errors it is desirable to have 
a uniform treatment for all the graphs, the number of which can be quite formidable. 
Secondly, the integrations over the Feynman parameters must be feasible, in the sense 
of fast convergence and high precision. 

At the two-loop level, we developed in ref. [^ a framework by which all two- 
loop processes with general renomalizable interactions can be written as a small set 
of Feynman parameter integrations over ten special functions hi. It would seem 
that an expansion around {sq} should be an easy task. However, this is not as 
straightforward as it appears, because the scalar integrals are made of functions of {s}, 
and a brute force expansion will result in rather unwieldy expressions. What is needed 
is a procedure to expand the momentum integrals around {sq} before the momentum 
integration, in such a way that the final product is still a linear combination of the 
ten scalar integrals. This will greatly facilitate systematic and rapidly convergent 
numerical work. We shall describe such a procedure in the following, for two-point 
functions. 

Our method is based on an expansion of individual propagators around a finite 
external momentum. In the following, we use a simple one-loop example in QED to 
show how the procedure is carried out. Afterwards we extend it to two- loop and apply 
it to the b quark self-energy, and derive the wave function renormalization constant 
as needed for the Z ^ bb process. 

Let us consider the self-energy of a fermion (the discussion for the bosonic case 
proceeds similarly): 

j:{k) = [A{e)+^,A,{e)]^-k + B{e)+^,B,{e) . (1) 

We want to obtain an expansion of this expression in powers of k"^ around a finite 
value k"^ = k^. Given a set of Feynman graphs of a certain loop order, in the case when 
one is able to evaluate them analytically, it is of course straightforward to perform a 



momentum expansion of the functions in the above expression. 

However, when several finite masses and momenta are involved, in general such 
analytical solutions are not available beyond one loop. Instead, we can obtain sys- 
tematically an expansion of the integrands by expanding the individual propagators 
which contain the external momentum k. We introduce a deviation A^ of the external 
momentum fc^ around the finite expansion point ko^: 

k^ = fcoM + A^ (2) 

such that the deviation is orthogonal to ko, A- ko = 0. We have the freedom to make 
this specific choice for A because the coefficients in the expansion are independent of 
the direction of differentiation. Squaring k^, one obtains A'^ = k"^ — k^. 

When writing down the integrand of a Feynman graph, the numerator is already 
a polynomial of the external momentum k, while the propagators containing k can 
be expanded around ko: 



A^ + 2p-A [A^ + 2p ■ A]^ 



{p + ky — m? {p + koY — w? [{p + koY — rn?Y [{p + ko)"^ — m 

where p is a loop integration momentum, and we used the fact that A ■ /cq = 0. 

Thus the dependence of the integrand on the deviation A becomes polynomial. 
For each order in A one can further use the algorithm described in ref. [|l| to calculate 
directly the desired coefficient. One obtains directly the expansion of the functions 
A(fc^), A^iji^), B{k'^) and B^{k'^). In particular, one is able to extract directly the 
mass and wave function renormalization counterterms. 

As a simple example, let us look at the electron mass operator at one-loop order 
in quantum electrodynamics. The momentum space integral in Feynman gauge is : 

,yil,^ = .2 /■ dy 7^[7-(P + fc) + ^]7A. ^,. 

in which m is the electron mass, and m^ is the photon mass regulator. The conven- 
tional way of introducing a Feynman parameter x and integrating over p gives: 

S(fc) = t/tit /rfxr(-i)(m2)f [nm + (2 - n)(l - x)j ■ k], (5) 

[Attj 2 J I 

where m^ = xm^ + (1 — x^w? — x(l — x)A;^, and n = 4 + e. This can be expanded in 
e, and after a MS subtraction we have: 



S(^) = -^—^ I dx\m\2x^1{\^x)\^{ml)\ 

^rr?x{\ — x^\ 



■{'-f ■ k — m) 



;i-x)[2 + 21nK)] + 



ml 



+0((7-fc-m)2)} , (6) 
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Figure 1: The two-loop Feynman graphs which contribute to the top-dependent b quark 
self-energy at 0{g'^as). 

with m"^ = x^m? + (1 — x)m'i. 

In the second approach we just proposed, we expand the electron propagator 
according to eq. 3, and write: 



S = So + El + S2 + 0(A=^ 



(7) 



where {p = p' — xko): 



iSo 



(2vr) 



dxd'^p' 



,„ ,nm + (2 - n)7 ■ [(1 - x)ko + A] 



(2vr) 



dxd'^p X 



(p'2 — miy 
, 2(^y27 ■ A + A2[nm + (2 - n){l - x)-f ■ k^ 



iSs = 87;^^ / dxx^dy'Ap'^A^ — 



(27r) 



(j)l2 _ 7722)3 

+ (^)(l-a:)7-fco 

(•p/2 _ 771^)4 



After making a MS subtraction in eqs. 8 and adding up the resuhs, we regain 
the result in eq. 6. This confirms the consistency of the expansion procedure. Note 
the terms 7 ■ A in eqs. 8, which combine with the 7 ■ k^ terms to produce 7 ■ k terms 
in eq. 6, as they must. This feature can be used as a check. 

This expansion procedure can be extended directly at two-loop level. The resulting 
coefficients of the A expansion can be further decomposed into hi functions and 
evaluated numerically along the lines of ref . [^ . 



As an example, we give in figure 1 the two-loop Feynman graphs which contribute 
to the top-dependent b quark self-energy at 0{g'^as)- The wave function renormal- 
ization constant derived from this two-point function contributes to the important 



process Z ^ bb [H-|10 



By expanding the propagators according to eq. 3, and identifying the terms 
according to their power of A, one can obtain the momentum expansion of the func- 
tions A{k'^), A5{k'^), Bik"^) and B^ik"^) of eq. 1 on-shell, i.e. around k^ = ml. The 
wave function renormalization contribution involves A{ml), A^lml), B'{ml), B'^{ml), 
A' {ml), and A'^{mLl), where the prime denotes the derivative with respect to the mo- 
mentum squared {d/dk"^). 

After isolating from the diagrams the A{ml), A^{ml), B'{ml), B'^{ml), /l'(m^), 
and A'^{mLl), contributions, we use the general algorithm of ref. [|I| to evaluate them. 
We used Mw = 80.41 GeV, m^ = 4.2 GeV, and sin^ 6'vk = .231. To regularize the 
infrared singularities, we introduced a mass regulator for the gluon. This being an 
0{g'^as) correction, the use of a gluon mass regulator is legitimate. When calculating a 
physical process, the gluon mass singularities are canceled by those from real emission 
diagrams; the net effect is to replace the mass regulator by the experimental energy 
resolution. The final results are given in figure 2 for a range of the top mass. Here 
we used MS to treat the ultraviolet infinities. 

At this point we note that in ref. p| the b quark self-energy is given in the massles 
b limit as a large top mass expansion up to the lO*'' order. A direct comparizon of 
these results with ours is difficult because ref. uses dimensional regularization for 
handling the infrared divergencies, while we have a gluon mass regulator. However, 
we checked that in the massles b limit our result for the self-energy contains only 
a left-handed component {i.e. A{k'^) = A5{k'^), B{k'^) = B^{k'^) — > in eq. 1), in 
agreement with the results of ref. ||^ . 

To conclude, we have shown how momentum expansions of massive two-loop dia- 
grams around a finite external momentum can be obtained systematically and eval- 
uated with the general methods which we introduced previously in ref. IQ. Such 
momentum expansions are needed in the computation of radiative corrections, for 
instance in the evaluation of wave function renormalization constants. The method 
we described in this letter has the advantage that provides a uniform treatment of all 
two-loop topologies, and is suitable for implementation in a computer algebra pro- 
gram for an automatic treatment. The final result of the computer algebra program 
is a set of special functions hi which are evaluated further numerically. 
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